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ABSTRACT 

We searched through roughly 12 years of archival survey data acquired by the Katzman 
Automatic Imaging Telescope (KAIT) as part of the Lick Observatory Supernova 
Search (LOSS) in order to detect or place limits on possible progenitor outbursts of 
Type Iln supernovae (SNe Iln). The KAIT database contains multiple pre-SN images 
for 5 SNe Iln (plus one ambiguous case of a SN Iln/imposter) within 50 Mpc. No 
progenitor outbursts are found using the false discovery rate (FDR) statistical method 
in any of our targets. Instead, we derive limiting magnitudes (LMs) at the locations of 
the SNe. These limiting magnitudes (typically reaching mji r; 19.5 mag) are compared 
to outbursts of SN 2009ip and r] Car, plus additional simulated outbursts. We find 
that the data for SN 1999el and SN 2003dv are of sufficient quality to rule out events 
40 days before the main peak caused by initially faint SNe from blue supergiant 
(BSG) precursor stars, as in the cases of SN 2009ip and SN 2010mc. These SNe Iln 
may thus have arisen from red supergiant progenitors, or they may have had a more 
rapid onset of circumstellar matter interaction. We also estimate the probability of 
detecting at least one outburst in our dataset to be > 60% for each type of the example 
outbursts, so the lack of any detections suggests that such outbursts are either typically 
less luminous (intrinsically or owing to dust) than ^ —13 mag, or not very common 
among SNe Iln within a few years prior to explosion. 

Key words: supernovae: Type Iln — luminous blue variables — supernova progen¬ 
itors 


1 INTRODUCTION 


Core-collapse supernovae (CCSNe) exhibit a wide diver¬ 
sity in their light curves and spectral properties, and it 
remains an enduring challenge to connect these properties 
to the evolution of their progenitors. Prior to 2009, four¬ 
teen CCSN progenitor stars or progenitor outbursts (eight 
SNe II-P, three Ilb, one Iln, one Ibn, a nd one Il-pec) had 
been detected in pre-explosion images (ISmartd l2009l . and 
references therein). Over the last few years, at least six 
more (Iln: 1961V, 2009ip, 201 011; lib: 2011dh, 2013df; Il n- 
P: 2011ht) have been found dSrn ith et all |2010al, |2011alf9; 
Kochanek et al.l 2011 ; Folev et al 20 111: Maund et al.l 20 111: 


Van Dvk et al.ll2014l : Fraser et al. 2^^J . Detection of these 

progenitors is crucial for gaining a better understanding 
of the observed SN diversity. Some theoretical models for 


* E-mail: cgbilinsk@gmail.com 
f Deceased 2011 December 12 


the end stages of massive-star evolution suggest that stars 
with 30-60 Mq will collapse to form a black hole with¬ 
out a SN explosion (iFryeil 1 19991 : iHeeer et al.l l2003l . but 
see lO’Connor &: Otu l2013l l. The detection of what seem 
to be massive luminous blue variable stars (LBVs) as 
the progenitors of the Type Iln SNe 1961V ISm ith et all 
2011^ Kochanek_et_^ 2ClN), 2005gl (|Gnb^ni_et_^ 2007; 


GnLYamfc LeonardI |2009 Il 2009ip (ISrnith et al.l |20~ : 
Folev et al. I I2OIII I. md 2010jl llSmith et al.ll2011lj ). however. 


challenge this notion. 

Optical spectra of SNe Iln exhibit prominent, relatively 
narrow hydrogen emission lines which arise from the inter¬ 
action of the expanding SN debris with a recently ejected 
circumstellar shell that is typica lly moving at < 1000 km s“^ 
llSchlege]||l990l : lFiliDpenk^ 1997fl . Because these lines appear 
soon after the SN explosion, it is expected that the shell of 
circumstellar matter (CSM) was ejected by t he star withi n 
a few years before the SN explosion itself dSmi thl l2014l'l . 
Such pre-SN eruptions in the years just before the SN sug- 
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2 Bilinski et. al. 


gest that instabilities arise dur i ng the late stages of nucl ear 
burning (ISmith fc Arnetjl2014l : [Quataert fc Shiod3l2012l i to 
power mass-loss episodes that might result in a detectable 
brightening of the progenitor system. In fact, there have 


SN 2006jc fibn: [Folev et al. 

2 OO 7 I: jPastorello et all 

2007h. 

SN 2009ip (Iln: ISmith et al 

2010a]: iPastorello et al. 

I 2 OI 3 I: 

Mauerhan et all l2013a]l, SN 2011ht (Iln-P: Fraser et all 

2013a]L and arEuablv SN 1961V flln: ISmith et al.ll2011bl: 

Kochanek et all I 2 OII]). thus 

strengthening the theoretical 
nd SNe Ilnjj Even when no 
ade, inferred wind and mass- 
. 20071: Gal-Yam et al. 20071: 

connection between LBVs 
direct detection has been m 
loss parameters llSmith et a 


20131: ISmithl l2014h. strong infrared excesses at late times 


i Gerardv et aLll2002l : lFox et al.ll2011^ . and a theoretical con- 
necti on to quasi-periodic radio modulations (iKotak fc VinkI 
I 2 OO 6 II su ggest that LBV-lik e progenitors are likely for some 
SNe Iln. lofek et al.l (l2014al l estimate that > 50% of SNe Iln 
in their Palomar Transient Factory (PTF) data have pre¬ 
explosion outbursts brighter than M = —14 mag, though 
some of these pre-peak brightening events may actually be 
the SN explosion itself (this is discussed more in M.4I) . 

Recent work suggests that SNe Iln follow LBV-like 
eruptions, although it remains unclear if the driving insta¬ 
bility of the pre-SN eruptions is t he same as t he instability 
of classical LB Vs (see review by ISmi thll2014l L In the con¬ 
text of stellar evolution, LBVs were suggested to be a tran¬ 
sitional phase between 0-type stars and Wolf-Rayet stars 
that skipped a red supergiant phase. Giant eruptions, oc¬ 
curring w hen a massive star inc r eases its bolometric lu¬ 
mino sity dHumphrevs et all 1 19991 : iHumphrevs fc DavidsorJ 
Il994h . were included in the LBV phenomenon. By spatially 
resolving the circumstellar shells resulting from two such 
LBV outbursts within our own Galaxy (P Cygni and r] 
Garinae), studies have been able to measure the amount 
of mass lost in observed LBV giant eruptions. Estimates 
of P Cygni’s circumstellar shell from the 1600 AD out - 
burst total only about 0.1 Mq dSmith fc Hartig^ l2006ll . 
whereas estimates of rj Carinae’s shell from its mid-19th 
century outburst total around 10-20 M m dSmith et al1l2003l : 
ISmith fc Ferlan3 l2007l : I Gomez et al.l 20101 1. While other 
Galactic LBVs h ave been observed wit h GSM shells sur¬ 
rounding them dSmith fc Owocl^ I2OO6I I. only these two 
events were actually observed during the eruptions that pro¬ 
duced their GSM shells. 

A number of recent developments have challenged the 
traditional view of LBVs, including evidence of shock- 
powe r ed even ts rather than super-Eddington winds dSmi thi 
I2OO8I . l2013al l. light-e cho spectra that do not match ex - 
pectations for winds dRest et al.ll20lj : IPrieto et al.ll2014h . 
and a wi dening range of initial masses that experience 
outbursts dPrietdr2 0o4 iPr ieto et al]l2008l : iThompson et al.l 
I2OO9I I. [Smith fc Tombleson ( 2014h show that LBVs tend to 
be more isolated than O-type stars and even Wolf-Rayet 


^ Note that we distinguish between a detection of a progenitor 
and a detection of variability that indicates a pre-SN outburst. It 
remains possible that some SNe with only one epoch of a detected 
progenitor may have in fact been detected during a relatively 
bright outburst phase. 


stars, challenging the idea that they are in transition be - 
tween these two phases. Instead. ISmith fc TomblesonI d2014ll 
propose that they may be the late evolutionary stage of 
mass gainers or mergers in binary systems. Finally, as noted 
above, the discovery of a number of LBVs exploding as CC- 
SNe without first experiencing a relatively long (> 0.5 My) 
Wolf-Rayet phase also challenges the traditional view of 
LBVs. 

In this paper, we further explore the connection between 
LBVs and SNe Iln by searching for additional progenitor 
outbursts preceding the SNe Iln. The data were acquired 
starting in 1998 as part of the Lick Observatory Super¬ 
nova Search (LOSS) with t he 0.76 m Katzman Au tomatic 
Imaging Telescope (KAIT: iFilippenko et al]|2001^ . KAIT 
began automatic operati on and discovered its first SN in 
1997 (iTreffers et al.l[l997ll . While the progra m was primar¬ 
ily fo c used on the det ermination of SN rates dLeaman et al.l 
I2OIII: iLi et al.l l2011blR> and SNe la for use as cosmolog¬ 
ical probes ( Ganeshahngam et al.l l20ld . I2OIII I , the KAIT 
database includes pre-SN images and SN photometry of 
many SN event s that can be utilised for other purposes 
(iLi et al.l l2011al L SN 1997bs, the first object that KAIT 
discovered with a SN designation, was interpreted as not 
being a SN at all, and became a prot otype for the class 
of “SN imposters” (IVan Dyk et al.ll200Cll j. A recent study by 
lAdams fc Kochanek ( 2015l L however, raises some additional 
questions about whether or not this was a terminal explo¬ 
sion. The discussion of SN Iln progenitors has intensified 
because of the recent discovery of erupting progenitor sys- 
tems like SN 20 09ip (jSmith et al.l l2010al : iFolev et al.l I2OIII : 
iPastorello et al.ll2013lj. 

Here we aim to constrain the frequency and luminosity 
of such outbursts. Information about the KAIT/LOSS pro¬ 
gram and the targets selected for our study is presented in 
m An analysis of our statistical approach to searching for 
progenitor outbursts and then constraining the brightness of 
these events with limiting magnitudes is given in 33] Section 
13] discusses the nature of the progenitors and our fractional 
coverage rates. We conclude in 31] with implications for fu¬ 
ture searches. 


2 OBSERVATIONS AND DATA ANALYSIS 

2.1 Data Acquisition, Reduction, and Processing 

Unfiltered photometry of over 850 SNe has been obtained 
since 1998 with KAIT (0.76 m) at Lick Observatory. Over 
time, the telescope has used three different CGDs, each with 
a response that when combined with the optical path effi¬ 
ciencies results in a sensitivity most similar to th at of the 
R band (iLi et al.ll2003l : iGaneshalingam et alJl201ol l . For this 
reason, we do all of our photometric analysis and compari¬ 
son in R when possible. With a scale of 0'.'8 pixel“^ and an 
effective chip size of 500 x 500 p ixels, KAIT’s total field of 
view is 6'. 7 x 6'. 7 (iLi et al.ll2003l) . Typical integration times 
for KAIT span the range of 16-40 s. During much of the time 
since 1998, KAIT has monitored roughly 15,000 galaxies at 
redshifts 2; < 0.05 in search of transients. Over 1200 galax¬ 
ies can be observed in a single long winter night (but fewer 
than 800 in a short summer night), resulting in a cadence of 
roughly 3-10 days for most targets (excluding the portion of 
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Constraints on SN Iln Progenitor Outbursts from LOSS 3 


the year that the target cannot be observed because it is in 
the daytime or bright twilight sky). These transient candi¬ 
dates go through a process of being automatically compared 
to template images and flagged, verified by human checkers, 
and then reobserved the following night before a SN is an¬ 
nounced llFilipp enko et al.ll200ll : lOaneshalineam et akllioill : 
iLeaman et aLll201ll l. 

Reduction and processing of the KAIT images wa s per¬ 
formed using techniques described bv iLi etliD (l2011al l. Bias 
corrections and flatfielding are automatically performed at 
the telescope for each target. The primary template image 
used for alignment and galaxy subtraction was obtained as 
part of the routine observation program on a photometric 
night with a longer exposure time than the typical images 
(usually 40 s as opposed to 16-30 s). A deep template with 
the SN present was generated using all images up to one year 
after the SN explosion that were of high quality and still 
contained the SN event as verified visually after galaxy sub¬ 
traction. This deep SN template was used to provide precise 
astrometry of the SN for the later analysis using artificial 
star injection. 

We obtained astrometric solutions for our primary tem- 
plate image by up loading our image to Astrometry.net 
dBarron et al.ll200a '). We then applied these astrometric so¬ 
lutions to the data images after they were aligned with the 
template image. The IRAF tasks phot, psf (on one of the 
brightest stars away from the galaxy or edge of the image), 
and allstar in the DAOPHOT package provided aper¬ 
ture and point-spread-function (PSF) photometry for all 
of our images. We calibrated the measured magnitudes to 
known magnitudes for any stars also present in the USNO- 
B databa se, which has rath er large uncertainties of about 
0.25 mag dMonet et al.ll200^ 1. We performed image subtrac¬ 
tion using the hotpants prograrrQ after which we could 
use the DAOPHOT package addstar task to inject arti¬ 
ficial stars. The phot task was used again to measure the 
photometry of the artificial stars and confirm that these ar¬ 
tificial stars in fact matched the known magnitudes that we 
inserted. We combined our images by using the imcombine 
task in IRAF in average mode when needed. The combined 
images were often cropped in order to obtain an intersection 
of our data that excluded regions which only contained data 
from a number of images less than the total in the stack. 


2.2 Type Iln Supernova Targets 

From the large sample of KAIT-observed SNe, we first se¬ 
lected all SNe Iln, candidate SNe Iln, and SNe II as classi¬ 
fied in the lAU Central Bureau of Astronomical Telegrams 
(CBAT) list of SNe. We then restricted this sample to tar¬ 
gets within 50 Mpc because KAIT images having a limiting 
apparent magnitude of mu « 19.5 result in an absolute mag¬ 
nitude of Mr « —14 at this distance (neglecting reddening). 
This is the peak luminosity of r] Carinae’s 19th-century erup¬ 
tion to which we compare the progenitor outburst limiting 
magnitudes, and it is a typical va lue for SN imposte rs (al¬ 
though they exhibit a wide range: ISmith et ani2011al l. The 
classification of the objects in the lAU CBAT list of SNe 


^ hotpants is available at http://www. astro . Washington. 
edu/users/becker/v2.0/hotpants.html 


is usually based on the first announcement published, but 
these are often revised upon further study. Accordingly, we 
checked the literature for each of our targets to verify that 
they were indeed SNe Iln. Overall, this selection process pro¬ 
vided us with five unambiguous SNe Iln and one ambiguous 
SN Iln/SN imposter (SN 2006am). 

Figure [1] shows stacked images of each SN field with 
the SN location labeled. Table [T] enumerates many of the 
important properties of each of the targets. In all but 
one case (SN 1999el), use of Hubble’s law was the only 
method available for distance determination for our tar¬ 
gets. We obtain these redshift-based distances from the 
NASA/IPAC Extr agalactic Databas fl which assumes Ho = 
73kms“^ Mpc“^ (iRiess et al.ll200a i and takes into account 
influences from the Virgo cluster, the Great Attractor, and 
the Shapley supercluster. Extinction along the line of sight 
to eac h SN’s host galaxy due to the Milky Way was ap¬ 
plied (ISchlaflv fc Finkbeinerll201ll i. Extinction produced in 
the host galaxies themselves is not known for any of our 
targets. Since these SNe Hn show signs of prior mass loss 
resulting in CSM, it is possible that the progenitors were 
enshrouded in dusty CSM that was destroyed by the SN. 
Our limiting magnitudes do not account for this possible 
additional extinction. 


2.2.1 SN 1999el 


SN 1999el was discovered near NGC 6951 as a part of 
the Beiji ng Astronomical Observatory Supernova Survey 


(BAOSS; ICao et al.l 


1999|) and peaked at mn = 14.5 mag 


{Mr = —18.4 mag; Di Carlo et al.l 1^00211 on 8 November 
1999 (UT dates are used throughout this paper). It was de¬ 
term ined to be a SN Hn based on BAO spectra dCao et al.l 
I 1999 II . Eollow-up observations in the near-infrared also con¬ 
firm its classification as a SN Hn (iD’Alessio et al.l I2001I : 
IPi Carlo et al.l I2OO2I : iDudlev fc Eischa 1999l l. The ultravi- 
olet/optical luminosity resulting from interaction between 
the SN shock and pre-existing circumstellar dust caused a 
blueward shift of the J , JT, and K light cu rves within 5- 
80 days after discovery (ID’Alessio et al.ll2001^ . Light echoes 
from hot CSM dust are also thought to play a role in the 
near-infrared light curves obser ved over the course o f 416 
days following maximum light llDi Carlo et al.lf2002l L These 
signs of interaction between the SN explosion and CSM are 
likely the result of mass-loss episodes prior to the final explo¬ 
sion. The distance to SN 1999el (26.3 Mpc) was determined 
by us ing the Type la SN 2000E in NGC 6951 ilPi Carlo et al.l 
I2OO2II . A Gal actic extinction value of Ar — 0.808 mag was 
adopted from ISchlaflv fc Einkbeineil ([20l3). 


2.2.2 SN 2003dv 

SN 2003dv was discovered near UGC 9638 as part of the Lick 
Observator y and Tenagra O bservatory Supernova Search 
(LOTOSS: Ixotak et al.ll2003l ~) and peaked at an unfiltered 
apparent magnitude of 15.4 {Mr = —17.5 mag) on 4 May 


^ The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech¬ 
nology, under contract with the National Aeronautics and Space 
Administration (http: //ned. ipac .caltech.edu). 
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RA (J2000) 


SN 2006am 
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SN 2010jl 
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04.00s 02.00s OlmOO.OOs 58.00s 13h00m56.00s 

RA (J2000) 


Figure 1. Stacked images for each of the SNe Iln showing the SN location. Each image contains the actual SN except that of SN 2011A 
because no images were taken with the SN present. For SN 2011A, the tick marks indicate where the SN would have been had we 
obtained post-SN images of this region. 


200 fl Spectra taken about a week before maximum bright¬ 
ness suggest that it was a SN Iln based on its narrow 
(350kms“^) Hq profile, which is slightly redward of the 
peak of th e broad component (10, 000 kms“^) of this Balmer 
emission (iKotak et alj [iool l. The distance to UGC 9638 
obtained from NED as described in H2.2I is 38.0 Mpc. A 

^ Light curve taken from http: / /www. rochesterastronomy. 
org/sn2003/sn2003dv.html 


Galactic extinction of An = 0.030 mag was adopted from 
ISchlaflv fc Finkbeineil (l201ll ~). 


2.2.3 SN 2006am 

SN 2006am was discovered near NGC 5630 as a part of 
LOSS with an unfiltered apparent magnitude of 1 8.5 {Mr = 
— 14.8 mag) on 22 February 2006 jLee fc Lill2006lf . This ab¬ 
solute magnitude indicates a rather low luminosity for a 
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Table 1. List of Type Iln Supernova Targets and their Properties 


SN 

Host Galaxy 

UT Peak Date 

RA (12000)=' 

Dec (J2000) 

Dist. (Mpc)^ 

Peak App. Mag. (Band) 

Peak Abs. Mag. 

1999el 

NGC 6951 

Nov-08-1999 

20'"37"* 17.72= 

-|-66°06'11.5" 

26.3‘= 

14.5 (R) 

-18.4 

2003dv 

UGC 9638 

May-04-2003 

14 58 04.92 

-1-58 52 49.9 

38.0 

15.4 (uf/R) 

-17.5 

2006am 

NGC 5630 

Feb-24-2006 

14 27 37.24 

+41 15 35.4 

44.2 

17.0 (uf/R) 

-16.2 

2008fq 

NGC 6907 

Sep-22-2008 

20 25 06.19 

-24 48 27.6 

47.4 

15.4 (uf/R) 

-18.1 

2010jl 

UGC 5189A 

Nov-05-2010 

09 42 53.33 

+09 29 41.8 

48.9'* 

12.9 (uf) 

-20.6 

2011A 

NGC 4902 

Jan-2-2011 

13 01 01.19 

-14 31 34.8 

38.2 

16.9 df) 

-16.1 


"^Peak dates are based off of published data if a well-sampled light curve exists. In the case that the published data are limited, we 
instead determine the peak date from the KAIT data. 


^Unless specifically cited, distances were obtained from the NED database, which assumes Hq = 73kms~^ Mpc“^ and takes into 
account influence from the Virgo cluster, the Great Attractor, and the Shapley supercluster. Milky Way extinction values along the line 
of sight to the host galaxies were taken into account when determining absolute magnitudes llSchlaflv Si Finkbeinedboilli . 


toi Carlo et al.l ll2002li . 
^Smith et al.l <l2011lji . 


CCSN, suggesting that SN 2006am may be either a SN 
imposter or a significantly extinguished SN. However, SN 
2006am exhibits properties of a SN Iln with a blue con¬ 
tinuum and Balmer emission lines containing both a nar¬ 
row component (300- 400 kms~^) and a weaker broad com¬ 
ponent (2000 kms~^: [ijee et al. lEooi). Given its faint mag¬ 
nitude along with its blue continuum and narrow spectral 
lines, SN 2006am may well be a SN imposter — but it may 
also be a SN Iln, so we include it in our sample. The dis¬ 
tance to NGC 5630 obtained from NED as described in ^2.21 
is 44.2 Mpc. A Galactic extinction of Au = 0.025 mag was 
adopted from ISchlaflv fc Finkbeineil ([^OlJ). 

2.2.4 SN 2008fq 

SN 2008fq was discovered near NGC 6907 as a part 
of LOSS with an unhltered apparent ma gnitude of 15.4 
(Mu = —18.1 mag) on 15 September 2008 (IThrasher et al.l 
l2008h . Spectra showed narrow and blueshifted (670kms“^) 
P-Cygni-like Ha and H/3 components superposed on 
a blueshifted (75 00 km s~^) diffuse absorption feature 
ll Quinn et al.l l2008l l . Like wise, the Na I line shows a sim¬ 
ilarly complicated prohle. lOuinn et al.l ll2008ll suggest the 
presence of multiple shells in the SN ejecta arising from 
ejecta-CSM interaction. Although originally classihed as a 
Type II SN, the narrow Ha features overlaying a broader 
component probably indicate that this object was a SN Hn. 
The distance to NGC 6907 obtained from NED as described 
in il2.2l is 47.4 Mpc. A Galactic extinction of Au = 0.137 mag 
was adopted from ISchlaflv fc Einkbeineil ll201 ih . 


2.2.5 SN 2010jl 

SN 2010jl was discovered near UGC 5189A by 
iNewton fc PuckettI ll2010lj and peaked at an unhltered 
apparent magnitude of 12.9 (M = —20.6 mag) on 5 Novem¬ 
ber 2010 (ISmith et al.l l2011bl L This absolute magnitude 
places SN 2010jl in the class of sup erluminous super novae 
(SLSNe ) along with SN 2003ma (iRest et al.l 2011 ), SN 
2006tf (ISmith et al.l l2008l j. and SN 2006gy l Ofek et ^ 
[20 o 3). SN 2010jrs host galaxy, UGC 5189A, is reported to 
be moderately metal poor with Z ~ 0.2-0.5 Zq (Istoll et al.l 


[2^. Spectra taken two weeks af ter discovery exhib it a 
narrow component of ~ 120hms”*^ (ISmith et al.ll2011lJ) . 

Archival Hubble Space Telescope (HST) data reveal a 
luminous blue point source at the loca tion of SN 2010jl ap - 
proximately 10 yr before its explosion (ISmith et al.ll2011bll . 
Since the SN has not yet faded, this progenitor candidate 
could be a massive young star cluster, a quiescent lumi- 
nous blue star, or a s tar undergo i ng an LBV-like eruption 
dSmith et al.l l2011blj . lOfek et"^ (l2014bll combined visible 
and X-ray data to estimate the presence of > 10 Mq of GSM 
surrounding the progenitor of SN 2010jl. Because of its ex¬ 
treme brightness and relative proximity, SN 2010jl has been 
the focus of numerous further s tudies containing additional 


optical and X-rav observations dRov et al. 201ll: Stoll et al. 

2011 

: iPatat et alJ I 2 OIII: 

Andrews et al.l 2011: ISmith et al. 

2012 

:IChandra et al.ll2012 

:IZhang et al.ll2012l: Zoglauer et al. 

2013 

: Fransson et al. 20131). The distance to UGC 5189A 


obtained from NED as described in is 48.9 M£c, con¬ 
sistent with the value adopted bv ISmith et al.l (l2011blj . A 
Galactic extinction of Au = 0.059 mag was adopted from 
ISchlaflv fc Einkbeined |20l3). 


2.2.6 SN 2011A 

SN 2011A was discovered near NGC 4902 as a part of the 
Chilean Automatic Supernova search (CHASE) project at 
an unhltered appare nt magnitude of 16.9 (M = —16.1 mag) 
on 2 January 2011 (IPignata et al.l[2011^ . Spectra were in¬ 
terpreted to s uggest that SN 2011 A was a SN Hn similar 
to SN 2005cl (iPignata et al.l[2011^ . The distance to NGC 
4902 obtained from NED as described in SI2.2l is 38.2 Mpc. 
A Galactic extinction of Au = 0.109 mag was adopted from 
ISchlahv fc Einkbeinerl |2oT3). 


3 ANALYSIS 

Our primary objective was to detect or place upper limits on 
the possible progenitor outbursts at the locations of the SNe 
in our sample, within the time frame monitored by KAIT 
prior to the SNe. No progenitor outbursts were detected to 
a statistically signihcant level using the false discovery rate 
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6 Bilinski et. al. 


(FDR) statistical method (discussed in mB Therefore, we 
also used a number of different techniques to place limiting 
magnitude constraints for these SN progenitor outbursts as 
described below. 


3.1 False Discovery Rate Statistical Method 

The FDR statistical method dBeniamini fc Hochberdll995l l 
was used to study the possibility of a source being detected 
in the archival KAIT images for each of the SN targets. It 
controls for the fact that with a very large dataset, Gaussian 
fluctuations would result in some significant false detections 
if a naive hypothesis testi ng method was used even when no 
emission is truly present (iLi et al.ll2011al l. In order to apply 
this method, we first compute a p-value for the flux level 
in each image at the location of the SN. We then place the 
p-values in ascending order and attempt to find the largest 
k such that Pk ^ —a, where k is the index of the p-value, 
m is the total number of images considered, and a is the 
significance requirement. We take the null hypothesis to be 
that no emission was present at the locations of the SNe in 
the archival images. All images with k = *,..., ^largest are 
declared as rejections of the null hypothesis, which means 
that they are considered detections of a source. We chose a 
relatively low significance level of q = 0.05, which is roughly 
the 2a detection threshold, in order to see if there were any 
weak detections we could explore further. 

With these techniques, we found no statistically signif¬ 
icant sources of emission in any of our targets. To verify 
that the statistical approach was working correctly, we also 
ran the FDR test on SN 2010jl images within two months 
after the SN explosion and were able to obtain significant 
detections in each image. However, since we are not able to 
directly detect any progenitor outbursts, we focus instead 
on setting limiting magnitudes at the SN position for each 
of the targets and compare them to the well-known light 
curves of SN 2009ip and rj Carinae. rj Carinae serves as an 
example of a nonterminal eruption over extended periods of 
time, whereas SN 2009ip serves as an example of a terminal 
event which had brief outbursts within the years prior to ex¬ 
plosion. Although the nature of SN 2009ip’s 2012a event was 
initially debated (IPrieto et al.l [20131: l^ker fc Kashil l2013l: 


Pastorello_^t_ah 2013: [M auerha n et al.l ' 2013a HPraser et al.l 

2013bl : iMargutti et alJ l2014l l. recent evidence in favor 


of the 2012a event being a terminal explosion has not 
been disputed sinc e (Smith, Mauerhan, fc Prieto 2014, 
iGraham et al.l l2014l : iMauerhan et al.l 120141 ) . While both rj 
Carinae and SN 2009ip have extensive light curves that we 
can reference, we do not claim that all SNe Iln have out¬ 
bursts similar to those of these two events. As such, we also 
compare the limiting magnitudes we set to an array of sim¬ 
ulated outbursts that cover a wide range of parameter space 
in varying magnitudes, durations of outbursts, and number 
of outbursts in a given event. 


® We also searched for outbursts in SN Iln imposter events such 
as SN 2001ac and SN 2006bv, but found no statistically significant 
detections in this sample either. 


3.2 Limiting Magnitudes on Individual Images 

3.2.1 Artificial Star Injection 

We take two approaches to setting limiting magnitudes for 
our images. Our first approach involves the use of artificial 
star injection. After subtraction of a template image not 
containing the SN, we implant increasingly fainter stars in 
steps of 0.1 mag at the location of the SN, starting from 
mu « 15 mag until we no longer detect a source present at 
the location of the SN to the 3a level (merr < 0.362 mag) of 
our known inserted star. 

While this method should be the most realistic at set¬ 
ting a relevant limiting magnitude for our images, there are 
a few drawbacks. In the case where the SN is in the bright¬ 
est regions of its host galaxy, as is the case for SN 2008fq 
and SN 2010jl, subtraction errors become more significant, 
resulting in a lower threshold of sensitivity and a large in¬ 
crease in the standard deviation of our limiting magnitudes. 
This technique is also susceptible to error from cosmic rays 
in the signal flux aperture. Cosmic rays in the signal flux 
aperture cause systematically fainter limiting magnitudes 
to be set, though they are relatively infrequent in the data. 
Cosmic rays present in the sky annulus used for noise consid¬ 
erations are controlled with a standard rejection procedure. 
Lastly, we increase the magnitude of the artificial stars in in¬ 
crements of 0.1, causing some artificial scatter in our results 
owing to this discretisation. Results from this approach are 
shown in Figures 2-7 as red triangles and given in Tables 
2-7 as “Artificial Star LM, Individual.” 


3.2.2 Sky Background Noise Calculations at the SN 
Location 

The second way in which we determine limiting magnitudes 
for KAIT images is by measuring the noise level in a bias-, 
dark-, and template-subtracted image in an annulus around 
the SN location and calculating the flux that would be re¬ 
quired to produce a magnitude error of 0.362. This allows 
us to solve for the limiting flux required to reach this signal- 
to-noise ratio (S/N) given the sky parameters near the SN, 
without actually using the flux at the SN location itself. 
We convert this limiting flux to a limiting magnitude and 
compare the result to the artihcial star injection method. 

Although this method is subject to galaxy-subtraction 
errors in the sky annulus, it is not susceptible to cosmic rays 
nor galaxy-subtraction artifacts in the signal flux aperture 
because it does not use the flux in this region in determining 
the limiting magnitude. Gonsequently, it tends to produce 
less scatter. Even though this approach does not use the 
flux information exactly at the location of the SN, it does 
use the same annulus for background noise information as 
the artificial star approach. Results from this approach are 
given in Tables 2-7 as “Background LM, Individual.” 


3.3 Limiting Magnitudes on Monthly Stacked 
Images 

Since we expect the signal from a progenitor outburst to be 
quite faint (i.e., Mr « —14 mag) , we also took the approach 
of averaging all of the images within the same month in or¬ 
der to obtain deeper stacked images. These images were first 
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registered to a template image in preparation for stacking. 
After combination, we then ran the same artificial star in¬ 
jections and noise calculations at the location of the SNe. 
The stacked images often result in deeper limiting magni¬ 
tudes being set, but their temporal coverage restricts their 
utility in searching for outbursts on the ~ 10 day timescale. 
The primary utility of these stacked images is in providing 
coverage of the faintest [Mr = —13 mag) outbursts that we 
simulate. Results from this approach are shown in Figures 
2-7 as blue squares and given in Tables 2-7 as “Artificial 
Star LM, Stacked” and “Background LM, Stacked.” 


3.4 Results 


As mentioned in mi the FDR statistical approach resulted 
in no significant progenitor outburst detections. Therefore, 
we focus in this section on li miting magnitudes and opt not 
to calculate precursor rates as lOfek et al.l (|2014al ') do because 
our data only allow us to place constraints on the properties 
of the possible outbursts. Plots showing the limiting magni¬ 
tudes for each of the targets are included in Fig ures 2-7. The 
absolute magnitude light curves of SN 2009ip dSmith et al.l 
l2010al : iPastorello et al.ll2013l : iMauerhan et al.ll2013al l and n 
Carinae ( Smith fc Frewl 201ll l are included for comparison 
in these figures. The light curves shown are shifted in the 
time axis so that their peak observed magnitudes occur at 
t = 0. All of the limiting magnitude values for each of the 
targets are included in Tables 2-7. 


3 . 4.3 SN 2006am 

Results for SN 2006am are shown in Figure |4] and Table H) 
SN 2006am’s dataset covers many years, but it contains only 
two images within the months prior to the SN explosion. 
Because of its relatively faint [Mr — —14.8 mag) discovery 
magn itude according to initial rough photometry IILee et al.l 
|200^ . SN 2006am may be very extinguished, or perhaps 
a SN imposter instead of a SN Iln. No light curve exists 
for SN 2006am other than our two KAIT measurements, 
which are the same images used earlier to discover this event. 
When performing more thorough photometry on these im¬ 
ages, however, we find a peak magnitude of Mr ~ —16, 
suggesting that SN 2006am is more likely a SN Iln than a 
SN imposter. Most of the limiting absolute magnitudes for 
SN 2006am lie in the range from —13.5 to —15. 


3 . 4.4 SN 2008fq 

Results for SN 2008fq are shown in Figure [5] and Table (5] 
Unfortunately, SN 2008fq is located near the brightest part 
of its host galaxy (NGC 6907) and at a relatively large dis¬ 
tance of 47.4 Mpc. This results in lower S/N and brighter 
limiting magnitudes for the SN 2008fq dataset. However, the 
large number of images acquired for NGC 6907 mean that 
we have very good temporal extent and coverage for this 
SN, enabling us to rule out precursor outbursts in the —14 
to —16 absolute magnitude range rather effectively for this 
SN in the years observed by KAIT. The light curve of SN 
2008fq shows a slower decline than that of SN 2009ip. 


3 . 4.1 SN 1999el 

Results for SN 1999el are shown in Figure [5] and Table [21 
While SN lOOOel’s dataset is relatively limited in duration 
because the SN exploded less than two years after KAIT be¬ 
gan operation, it does contain many limits just before peak 
brightness and is the closest SN Iln that we consider at 
26.3 Mpc. Consequently, the majority of the limiting mag¬ 
nitudes set for the SN 1999el data lie in the absolute mag- 
nitude range from —1 3 to —14 mag. SN 1999ers light curve 
dPi Carlo et al.ir2002l l matches that of SN 2009ip (from the 
2012 b peak onward) very well both in magnitude and in 
shape. SN 1999el was not observed as early in its rise to 
peak, but its post-peak light curve shows evolution similar 
to that of SN 2009ip. As discussed below, we can rule out 
an initially faint SN if it were similar to SN 2009ip’s 2012a 
event. 


3 . 4.2 SN 2003dv 

Results for SN 2003dv are shown in Figure[3]and TableO SN 
2003dv’s dataset includes five years of limiting magnitudes, 
but only a handful of good images were acquired for this 
target each year. Most of the limiting magnitudes fall in the 
range of —13 to —14 in absolute magnitude for SN 2003dv. 
SN 2003dv’s light curve reaches a similar peak luminosity 
but shows a much slower decline compared to that of SN 
2009ip. 


3 . 4.5 SN 2 OIO 3 I 

Results for SN 2010jl are shown in Figure [6] and Table (6] 
SN 2010jl is located in a bright part of its host galaxy 
(UGC 5189A), and it is also the farthest SN Iln we con¬ 
sider [d = 48.9 Mpc). The number of good KAIT images of 
this object acquired each year is small, making this dataset 
noisy and sparse, even though it does cover twelve years. 
With most of the limiting absolute magnitudes falling in 
the range of —15 to —16 mag, and none being set in the 
months prior to SN 2010jrs peak brightness, this dataset 
does not allow us to place strong constraints on the na¬ 
ture of SN 2010jrs progenitor. A candidate progenitor with 
M = —12.0 mag was, however, detected for SN 2010jl in 
the F300W filter on j7S' r/WFP C2 roughly 10 yr before peak 
brightne ss (ISmit h et al.ll2011bl b The light curve of SN 2010jl 
(IZhang et al.ll2012l ') shows a significantly higher peak bright¬ 
ness than that of SN 2009ip and also a much slower decline 
from peak. 

3 . 4.6 SN 2011A 

Results for SN 2011A are shown in Figure [7] and Table [7| 
The dataset for SN 2011A covers a large temporal extent 
and contains enough images each year to provide good cover¬ 
age. Unfortunately, no KAIT data are available in the 1.5 yr 
immediately preceding the SN explosion because of a disk 
failure resulting in the loss of data between September 2009 
and September 2011. Most of the limiting absolute magni¬ 
tudes for SN 2011A lie in the range —13 to —14. No light 
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Table 2. Limiting Magnitudes (LM) for SN 1999eP 


Year 

Month 

Day 

Artificial Star LM, Individual 

Artificial Star LM, Stacked 

Background LM, Individual 

Background LM, Stacked 

1998 

10 

24 

19.2 

19.2 


19.4 

19.8 

1998 

10 

28 

19.1 

19.2 


19.6 

19.8 

1998 

10 

31 

18.8 

19.2 


19.7 

19.8 

1998 

11 

2 

18.5 

19.9 


19.6 

19.8 

1998 

11 

5 

18.8 

19.9 


18.8 

19.8 

""The entirety of this table is available electronically. A portion is displayed here for reference. 



Table 3. Limiting Magnitudes (LM) for SN 2003dv“ 





Year 

Month 

Day 

Artificial Star LM, Individual 

Artificial Star LM, Stacked 

Background LM, Individual 

Background LM, Stacked 

1999 

1 

28 

18.7 

18.7 


19.7 

19.7 

1999 

2 

1 

18.6 

19.2 


18.9 

19.4 

1999 

2 

27 

19.8 

19.2 


19.8 

19.4 

1999 

4 

10 

19.7 

19.9 


20.0 

20.2 

1999 

4 

15 

20.0 

19.9 


19.9 

20.2 


“The entirety of this table is available electronically. A portion is displayed here for reference. 


Table 4. Limiting Magnitudes (LM) for SN 2006am“ 


Year 

Month 

Day 

Artificial Star LM, Individual 

Artificial Star LM, Stacked 

Background LM, Individual 

Background LM, Stacked 

1999 

1 

10 

19.5 

19.2 

19.3 

19.6 

1999 

1 

28 

18.5 

19.2 

19.6 

19.6 

1999 

2 

2 

18.3 

19.1 

19.5 

20.0 

1999 

2 

5 

19.6 

19.1 

19.4 

20.0 

1999 

2 

27 

17.2 

19.1 

19.5 

20.0 


•^The entirety of this table is available electronically. A portion is displayed here for reference. 



Table 5. Limiting Magnitudes (LM) for SN 2008fq“ 





Year 

Month 

Day 

Artificial Star LM, Individual 

Artificial Star LM, Stacked 

Background LM, Individual 

Background LM, Stacked 

2000 

6 

14 

19.6 

18.6 


19.5 

19.8 

2000 

6 

23 

17.8 

18.6 


19.4 

19.8 

2000 

6 

30 

18.2 

18.6 


19.6 

19.8 

2000 

7 

8 

19.0 

19.6 


19.5 

19.7 

2000 

7 

14 

19.1 

19.6 


19.0 

19.7 

•^The entirety of this table is available electronically. A portion is displayed here for reference. 



Table 6. Limiting Magnitudes (LM) for SN 2010jl“ 





Year 

Month 

Day 

Artificial Star LM, Individual 

Artificial Star LM, Stacked 

Background LM, Individual 

Background LM, Stacked 

1998 

11 

15 

18.0 

18.1 


19.4 

19.8 

1998 

11 

21 

18.0 

18.1 


19.6 

19.8 

1998 

12 

23 

17.2 

17.3 


19.8 

20.0 

1998 

12 

27 

17.2 

17.3 


19.5 

20.0 

1999 

1 

8 

18.0 

18.6 


19.6 

19.9 


“The entirety of this table is available electronically. A portion is displayed here for reference. 
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Figure 2. Archival KAIT limiting magnitudes for SN 1999el. Red triangular data signify limits set by data from single images, whereas 
blue-square data signify limits set by monthly stacks. (If only one image is available in a given month, the stack may consist of a single 
image resulting in overlap between the individual and combined limit.) The size of the combined data squares indicates how many images 
went into producing the combined image (the largest stack consists of 16 images and belongs to the SN 2008fq dataset). SN 2009ip’s 
light curve is overplotted as a dash-dotted green line, r] Carinae’s light curve is overplotted as a dashed black line. The vertical dashed 
black line marks the date of the peak observed magnitude. Near-peak light-curve values for SN 1999el were obtained from lPi Carlo et al.l 
1 I 2 OO 2 I I. We determine the date of the peak observed magnitude from this well-sampled light curve. 


Table 7. Limiting Magnitudes (LM) for SN 2011 A® 


Year 

Month 

Day 

Artificial Star LM, Individual 

Artificial Star LM, Stacked 

Background LM, Individual 

Background LM, Stacked 

1998 

12 

27 

19.1 

19.1 

19.0 

19.0 

1999 

1 

3 

16.9 

18.2 

19.3 

19.6 

1999 

1 

6 

18.2 

18.2 

19.2 

19.6 

1999 

1 

11 

18.4 

18.2 

18.7 

19.6 

1999 

1 

28 

17.7 

18.2 

18.1 

19.6 


“The entirety of this table is available electronically. A portion is displayed here for reference. 


curve is available for SN 2011A to be used in comparison to 
SN 2009ip’s light curve. 


4 DISCUSSION 

In general, the limiting magnitudes that we set prior to our 
sample of SNe Iln are not faint enough to rule out SN 2009ip 
or Tj Carinae-like eruptions with high confidence. Even the 
stacked images are usnally not snfficiently deep because 
the KAIT survey does not typically acqnire a large num¬ 
ber (> 10) of images of each particular galaxy each month. 
Thus, a main conclusion is that the lack of detections of pre- 
SN outbursts is not very surprising, even if all SNe Iln are 
core-collapse events with substantial {Mr « —14 mag) pre¬ 
cursor eruptions. Nevertheless, we can perform a number 
of statistical tests, which we discuss below, to gain a bet¬ 
ter understanding of what our limits imply with regard to 


SN Iln progenitor properties, and to guide future observing 
strategies. 

4.1 Coverage Rates for Example Outbursts 

Although control times were not calculated for the KAIT 
data for any SN imposters, we compute fractional coverage 
rates for hve different known LBV-like outbursts, which are 
shown in Figure [S] Three of these outbursts are taken from 
the SN 2009ip light curve (the last of which is actually likely 
to be the start of the SN e xplosion itself; see lMauerhan et al.l 
l2013al : ISmith et al.ll201^ ) and th e other two are taken from 
rj Carinae’s historical light curve (ISmith &: Frewll201 j ). The 
length of the outbursts varies from 45 to 192 days and some 
of the outbursts, like that of the second SN 2009ip-like burst, 
include many brief individual outbursts within them. Over¬ 
all, the outbursts cover a range from ~ —11 mag during a 
quiescent part of the burst to ~ —15 mag at peak outburst. 


© 2015 RAS, MNRAS 000.[Tlf22l 






























10 Bilinski et. al. 


All SN 2003dv Data 



12 

13 

14 

ic CD 

15 - 0 

16 - E 
O) 


17 


CD 


18 c 

Q 

19 S 
Q. 
Q. 

20 < 
21 
22 


Zoomed SN 2003dv Data 



Days 


Figure 3. Same as Fig. [2] but for SN 2003dv. The top plot shows all of the archival data while the bottom plot shows only data near 
peak magnitude. Near-peak light-curve values for SN 2003dv were obtained from the amateur astronomy data available on the Rochester 
Astronomy website of nearby SNe (http: //www. rochesterastronomy .org/sn2003/sn2003dv. html). We determine the date of the 
peak observed magnitude from this well-sampled light curve. 


consistent with th e observed distribu tion of SN imposter 
peak luminosities llSmith et al.ll2011a^ . Observations of SN 
2009ip were taken in either an unfiltered bandpass or the 
R band, while observations of rj Carinae were made primar¬ 
ily by eye in the visible. With the temperatures expected 
for these outbursts near peak 6000-7000 K), i?-band and 
visible photometry should be comparable (because the bolo- 
metric correction is small), except in the case of very high 
Ha emission levels. 

The fractional coverage rates we compute are essentially 
the probabilities that we would have been able to detect 


each particular outburst. We assume each of the outbursts 
has an equal likelihood of having occurred on any day since 
our first available observation for each SN. If the outburst 
is ever brighter than a limiting magnitude we have set for 
a given day, then we claim that we would have been able 
to detect that outburst. After inserting each outburst on 
every possible day that it could have occurred, we divide 
the total number of detected outbursts by the total number 
of outbursts inserted to find the fractional coverage rate. 
The results of this assessment are presented in Table [S] 
Because of the brightness of the host galaxies at the 
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Figure 4. Same as Fig. |3] but for SN 2006am. The two near-discovery light curve values for SN 2006am were obtained from our 
KAIT/LOSS survey data. We determine the date of the peak observed magnitude from this poorly sampled light curve. 


location of SN 2010jl and SN 2008fq, the limiting magni¬ 
tudes are too bright to place a strong constraint on the 
potential presence of pre-SN outbursts. However, the other 
four SNe have reasonable 0.10) fractional coverage rates 
for most of the example outbursts, with particularly high 
0.35) coverage rates for SN 2009ip-like burst ^3 simu¬ 
lations. Considering that we have no limiting magnitudes 
for a large portion of the year owing to right-ascension con¬ 
straints, any coverage rate near 0.5 is high. 

Even though we have the highest fractional coverage 
rates for SN 2009ip-like burst #3, we focus more on the 
other outbursts in this discussion because this event (2012a) 


is likely to be the SN explosion itsel f, which is initially faint 
because of having a BSG progenitor (iMauerhan et al]l2013al : 
ISmith et al.ll2of4 i. as was the case for SN 1987A. Our next 
most prominent outburst is that of rj Car-like burst #2 
(1843), which shows rates above 17% for all of our objects, 
excluding SN 2010jl and SN 2008fq. If all of the SNe Iln in 
our sample had gone through rj Car-like burst ^2 phases 
within their recent histories, then we would have an 88% 
chance of detecting at least one outburst from the six tar¬ 
gets. Even for SN 2009ip-like burst #2, the faintest and 
shortest duration outburst, we find a 64% chance of having 
detected at least one outburst among all of our events. 
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Figure 5. Same as Fig. [S] but for SN 2008fq. Near-peak light-curve values for SN 2008fq were obtained from our KAIT/LOSS survey 
data. Since only two data points exist on the light curve available from the Rochester Astronomy website of nearby SNe, we determine 
the date of the peak observed magnitude from our well-sampled KAIT light curve. 


4.2 Possible Extinction from Circumstellar Dust 

These calculations are heavily dependent on the fact that 
we ignore local extinction by circumstellar dust (though SN 
2009ip and r] Carinae may have been surrounded by dust 
as well). We consider the possibility of circumstellar dust 
enshrouding our targets by adding the effects of artihcial 
extinction to the example outbursts and recalculating the 
coverage rates. We find that with just Ar ^ 0.5 mag, the 
probability that we would have detected at least one out¬ 
burst from any of our targets similar to that of SN 2009ip’s 
burst ^2 drops below 10%. Thus, it may be that the ma¬ 


jority of SN Iln precursors are somewhat less luminous in¬ 
trinsically than SN 2009ip or rj Carinae, or that they appear 
fainter because of a modest amount of dust that formed in a 
previous eruption. The dust formed by precursor eruptions 
may be destroyed by the SN explosion. Because of this, we 
cannot use measurements of reddening along the line of sight 
to SNe to rule out extinction effects on the progenitor sys¬ 
tems due to CSM dust. 
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Figure 6. Same as Fig. [S] but for SN 2010jl. Near-peak light-curve values for SN 2010jl were obtained from l^ang et ^ lOd^i, 
determine the date of the peak observed magnitude from a light curve available on the Rochester Astronomy website of nearby 
which extends just a few days before the published light curve. 


We 

SNe 


4.3 Coverage Rates for Simulated Outbursts 

In order to further constrain the types of outbursts, we com¬ 
pute fractional coverage rates (the chance that we would 
have detected a giveu outburst if it had occurred ou auy 
random day since the start of our observations for the SN 
in question) for simulated outbursts with a variety of pa¬ 
rameters, for limiting magnitudes set by both the artificial 
star and background noise techniques for all of our SN Iln 
targets. The results of these tests are shown in Figures[9]and 
HqI We simulate outbursts similar to those from SN 2009ip 
(many short outbursts) and rj Carinae (less frequent but ex¬ 


tended outbursts) over a range of magnitudes. Specifically, 
in the case of SN 2009ip-like outbursts, we simulate ten-day 
light curves which vary in peak magnitude between —13 and 
— 16, and recur between one and ten times in a particular 
simulation. In the case of many repeated outbursts, the sep¬ 
aration between events is set to ten days. For instance, a 
five-event SN 2009ip-like light curve at m = —15 mag uses 
five separate outburst events that are ten days long at a con¬ 
stant magnitude of —15 with ten days separating each of the 
five outburst events. In the case of rj Carinae-like outbursts, 
we simulate a single light curve which varies in magnitude 
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Figure 7. Same as Fig. [S] but for SN 2011A. No plot for near-peak values is presented because no survey data for SN 2011A were 
available for these dates. No light curve was available for SN 2011A because of a disk failure resulting in the loss of data. 


Table 8. Fractional Probability of Detecting Given Outbursts Using Artificial Star Injection LMs 


SN 

Image Type 

2009ip Burst 

2009ip Burst #2 

2009ip Burst #3 

rj Car Burst ^1 

7] Car Burst ^^2 

1999el 

Individual 

0.48 

0.54 

0.56 

0.53 

0.62 

- 

Stacked 

0.20 

0.17 

0.46 

0.34 

0.40 

2003dv 

Individual 

0.32 

0.18 

0.48 

0.31 

0.47 

- 

Stacked 

0.23 

0.12 

0.43 

0.25 

0.45 

2006am 

Individual 

0.11 

0.01 

0.49 

0.01 

0.17 

- 

Stacked 

0.13 

0.02 

0.43 

0.03 

0.27 

2008fq 

Individual 

0.02 

0.00 

0.35 

0.00 

0.01 

- 

Stacked 

0.02 

0.00 

0.27 

0.00 

0.04 

2010jl 

Individual 

0.01 

0.00 

0.11 

0.00 

0.02 

- 

Stacked 

0.00 

0.00 

0.07 

0.00 

0.01 

2011A 

Individual 

0.15 

0.04 

0.40 

0.08 

0.26 

- 

Stacked 

0.14 

0.06 

0.36 

0.16 

0.26 

Collective® 

Individual 

0.74 

0.64 

0.96 

0.70 

0.88 

- 

Stacked 

0.56 

0.33 

0.92 

0.60 

0.83 


“These are the probabilities that we would have detected at least one outburst from any of the SN Iln targets. 


between —13 and —16, and duration of the outburst between 
ten and one hundred days. Note here that the SN 2008fq 
and SN 2010jl fractional coverage rates are diminished when 
the artificial star limiting magnitudes are used because of 
high host-galaxy brightness in the signal flux aperture. The 
remainder of the fractional coverage rates are quite high 
(> 0.2) for almost all cases plotted with M < —13.5 mag. 

The KAIT cadence of 3-10 days is sufficient for the 
purposes of this project, as can be seen from the varying 
duration of outburst plots in Figure [TOl As the duration of 
the outbursts decreases, the fractional coverage rates drop 
slowly compared to the steep dropoffs seen in the other plots. 
This slow drop is primarily caused by the fact that our lim¬ 
iting magnitudes only cover the ~ 50% of the year during 


which we have data, so that longer outbursts spend less time 
in our “blind spots.” 


Our main point is that the fractional coverage rates be¬ 
come very small for Mji > —13 mag, indicating that KAIT 
survey data are not sensitive enough to place constraints on 
outbursts fainter than this level. Because the SN 2009ip and 
T/ Carinae light curves spend most of their quiescent time 
at Mu « —11 mag or Mv ~ —12mag (respectively), such 
sensitivity levels are necessary to thoroughly constrain the 
frequency of their outbursts. 
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Figure 8. The top plot shows the entire light curves for SN 2009ip and r] Carinae, while the bottom is zoomed in to cover only ~ 3000 
days. The highlighted regions show the selected outbursts used for comparison to the calculated limiting magnitudes of each SN Iln. SN 
2009ip-like Burst #3 is actually the SN’s initial rise rather than a precursor outburst. 


4.4 Red vs. Blue Supergiant? 

Although it had been speculated for some time that red su¬ 
pergiants (RSGs) were the most likely progenitors of SNe II, 
evidence has suggested that some SNe II may explode while 
in a BSG or LBV phase. The most fa mous of such BSG 
progenitors is SN 1987A llArnettI Il989li . but more recent 
cases associated with SNe I ln hav e been suggested as well 
llSmith et al.l l2010al . l2011bl . I20l4l . The key difference be¬ 
tween a RSG and a BSG progenitor is the larger photo- 
spheric radius for a RSG progenitor. Because exploded stars 
with larger initial radii lose less thermal energy to adiabatic 
expansi on, their SNe can achieve a higher initial peak lumi¬ 
nosity (iDoggett fc Brandilll985f) . Whereas the RSG’s light 
curve essentially plateaus from this peak, a BSG’s light curve 
will continue to rise from its initially faint state as ^^ Ni decay 
begins to d ominate, as in the example of SN 1987A llArnettI 

Hiii,!!!!!). 

Interaction with CSM can also greatly increase the lu¬ 


minosity from a BSG (or any) explosion depending on the 
delay between the e xplosion time and the onset of strong 
GSM interaction (see lSmith et al.ll2014f) . but in this section 
we are concerned with the rise to peak brightness of the SN- 
ejecta photosphere. Since LBV progenitors to SNe Iln are 
expected to have BSG-like radii, we can look for these signa¬ 
tures in the light curve to constrain progenitor properties. 
SN 2009ip and SN 2010mc are strong candidates for hav¬ 
ing relatively compact BSGs prior to explosion (compact, at 
least, compared to a RSG of a similar luminosity). Because 
an abundance of data exists for both the SN explosion itself 
and the many outbursts in the preceding years, SN 2009ip’s 
light curve is used as a general template to explore our abil¬ 
ity to constrain the progenitors for the SNe Iln studied in 
this work. Given that some of our SNe Iln have limiting 
magnitudes set shortly before they reached peak brightness, 
we can use these limits to constrain how quickly the light 
curve rose in each case. For brevity, we refer to an initially 
faint BSG SN light curve as an 09ip/2012a-like event. 
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Coverage of one 10-day outburst with sky background LMs 


Coverage of one 10-day outburst with artificial star LMs 



Coverage of ten 10-day outbursts with sky background LMs 




Coverage of ten 10-day outbursts with artificial star LMs 



Coverage of /W=-14.0mag outbursts with sky background LMs Coverage of /W=-14.0mag outbursts with artificial star LMs 




Figure 9. Fractional coverage of SN 2009ip-like (many quick) simulated outbursts based on the limiting magnitudes for the various 
SNe Iln. The left column uses the limiting magnitudes based on background noise, whereas the right column uses the limiting magnitudes 
based on artificial star injection. The top row shows the fractional coverage for varying magnitudes of a 10-day simulated outburst. The 
middle row provides the fractional coverage for varying magnitudes of ten 10-day simulated outbursts. The bottom row gives the fractional 
coverage for varying numbers of 10-day outbursts at a fixed magnitude of M = —14. 
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Coverage of one 10-day outburst with sky background LMs Coverage of one 10-day outburst with artificial star LMs 



Coverage of one 100-day outburst with sky background LMs 


Coverage of one 100-day outburst with artificial star LMs 




Absolute Magnitude 


Coverage of one /M=-14.0mag outburst with sky background LMs 



Coverage of one M=-14.0mag outburst with artificial star LMs 



Figure 10. Same as Fig. (9] except for tj Carinae-like simulated outbursts. The middle row shows the fractional coverage for varying 
magnitudes of a single 100-day simulated outburst. The bottom row gives the fractional coverage for varying durations of a single outburst 
at a fixed magnitude of M = —14. 


SN 1999el, SN 2003dv, SN 2006am, and SN 2008fq have 
relevant data for this question. Figure [TT] shows zoomed-in 
light curves for these SNe near their times of peak bright¬ 
ness. With data ending just 14 days before peak, SN 1999el 
shows strong signs that it did not have a faint 09ip/2012a- 


like event prior to reaching its peak brightness. Data for SN 
2003dv are earlier relative to the peak, but do span an ef¬ 
fective range of 36 to 75 days prior to the peak brightness, 
suggesting a possible quick rise to peak as well. One limit is 
set 20 days before peak for SN 2003dv, but the noise level is 
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18 Bilinski et. al. 


too high to conclusively rule out a 09ip/2012a-like event. SN 
2006am also has limits 25 and 39 days before its observed 
peak (because we have so few data for the SN 2006am light 
curve, we do not know for certain that we actually have 
near-peak observations for this event), but the 25-day limit 
is from a poor-quality image and does not significantly con¬ 
strain the SN brightness. The observation 39 days prior to 
peak brightness does suggest that no detectable source at 
M « —14 mag was present at that time, but this alone does 
not rule out a 09ip/2012a-like event for SN 2006am. Pre-SN 
data near peak for SN 2008fq span a range of 16-135 days 
prior to the SN explosion, three of which seem to suggest 
a light curve fainter than that of SN 2009ip as it rose to 
peak. However, the majority of the limiting magnitudes are 
too bright to stringently constrain the rise to peak for SN 
2008fq. This is a direct consequence of its being far away at 
47.4 Mpc and being located near the brightest part of its 
host galaxy. 

Because SN 1999el and SN 2003dv do not seem to have 
had a faint pre-peak SN bump, it is unlikely that their pro¬ 
genitors were BSGs akin to SN 2009ip and SN 2010mc. In¬ 
stead, these SNe underwent a quick rise to peak that would 
be expected from a RSG progenitor, or perhaps from a BSG 
with immediate onset of strong GSM interaction. For SN 
2010jl and SN 2011A, we cannot place any constraints in this 
manner because no data are available within 150 days prior 
to peak brightness, but we do know that HS T images for 
SN 20 10]! suggest a possible BSG progenitor JSmith et al.l 

l2011bl) . 

We must also consider SNe Iln-P and their progenitors 
here because a SN Iln-P may appear as a SN Iln if no data 
are available to rule out a plateau phase. Hence, we can¬ 
not exclude the possibility that some of our SNe Hn are in 
fact SNe Hn-P. This is important because SNe Hn-P are 
hypothesised to arise from super asymptotic giant branch 
(AGB) stars which explode as electron-capture SNe, as in 


the case of SN 19 94W, SN 2009kn, SN 2011ht, an d the 
Grab nebula’s SN llChug ai et al.l 200'd: Kankare et al. 20ld; 


iMauerhan et al.l l2013bl : iFraser et al.ll2013al : ISmithll2013bll . 
SN 2011ht was reported to have an outburst o ne year prior 
to its terminal explosion llFraser et al.l l2013al l , suggesting 
that super-AGB stars may also produce nonterminal erup¬ 
tive mass-loss events just before exploding. The population 
of events that appear as SNe Hn may therefore have a diver¬ 
sity of progenitors from RSGs to BSG s rather than a single 
progenitor system llSmith et al.l 20091). 

Recent work by Ofek et al. ( 2014ah . which looked at 
a sample of sixteen SNe Hn in Palomar Transient Factory 
(PTF) data, claims that more than 50% of SNe Hn have 
at least one pre-explosion outburst brighter than 3 x 10^ L© 
(absolute magnitude M « —14) that occurs within 4 months 
before the SN0 This statistical estimate includes two events 
(SN 2010mc precursor 20 days before, and the PTF 12cxj- 
A precursor 10 days before) that occur just prior to the 


° This particular statistic assumes a homogeneous population 
among all of the SN Iln targets so that it can count multiple 
precursor events for a single target as essentially different events. 
Thus, if sixteen precursor outbursts were detected for a single one 
of the sixteen different targets, the statistic would claim that es¬ 
sentially 100% have such an outburst within the measured time 
frame. 


SN peak brightness, but are very similar to the 09ip/2012a 
prec ursor that occurred just 25 days before peak bright¬ 
ness. ISmith et al.l ll2013l l showed that the light curves of SN 
2010mc and SN 2009ip are almost identical. These “precur¬ 
sors” may well be the initial stages of the slowly rising SN ex¬ 
plosion itself be fore GSM interaction peaks, as was the cas e 
for SN 2009ip (IMauerhan et al.ll2013al : [Smith et al.ll2014l l. 
For these reasons, it seems that detectable luminous erup¬ 
tions (M < —14 mag) may be less common among SN Iln 
progenitors just be fore explosion than the 50% estimate by 
lOfek etlil] ( 2014al '). Our constraints on SN 2009ip-like and 
rj Carinae-like eruptions suggest that these more luminous 
eruptive events are uncommon, or that they are fainter than 
our example outbursts (intrinsically or owing to possible 
dust extinction). In fact, we find that our data require that 
< 40% of SNe Hn have pre-SN LBV-like eruptions of —13 
mag or brighter in the time frame covered by KAIT. 


5 CONCLUSIONS 

Our primary goal was to search through the KAIT archival 
database in order to constrain the properties of SN Hn pro¬ 
genitor outbursts by determining if any SNe Hn exhibit 
detectable outburst characteristics similar to those of SN 
2009ip or rj Garinae. No such outbursts were statistically 
detected in the dataset, implying that they are either typi¬ 
cally less luminous than expected or are not common among 
SNe Hn. Instead, we have provided limiting magnitudes for 
six SN Hn progenitors. 

With these limiting magnitudes, we are able to place 
constraints on the nature of a few of the progenitors to 
our SNe Hn. In particular, we find that SN 1999el and SN 
2003dv have constraining limiting magnitudes set within 40 
days prior to their peak brightness, suggesting that they are 
likely RSGs with steeply rising light curves or that GSM 
interaction turned on much more quickly than the delayed 
interaction in SN 2009ip and SN 2010mc. We also find that 
if the SN 2009ip and rj Garinae outbursts were typical of 
all SNe Hn, then we would have a good chance (> 60%) of 
detecting each given type of outburst at least once in our 
dataset. However, this is for an emerged luminosity corre¬ 
sponding to ~ —13 mag. This estimate ignores the possibil¬ 
ity that these objects (as well as SN 2009ip and rj Garinae) 
are in fact obscured by their own GSM dust, which is likely 
to be created if there have been previous such mass-loss 
episodes. It also depends somewhat on our particular view- 
ing angle for each SN if the GSM geometry is aspherical 
(IMauerhan et af]l2014h . 

In order to better study the nature of SNe Hn and their 
progenitors, it would be beneficial to perform surveys with 
deeper limiting magnitudes. The KAIT cadence of 3-10 days 
provides excellent temporal coverage of the outbursts we 
are exploring, while the (up to) 12 yr baseline of archival 
data provides nearly unparalleled temporal extent for such 
a project. However, the sensitivity of KAIT only allows us 
to reach limiting magnitudes of m « 19.5, which restricts us 
to SNe Hn within 50 Mpc. Any objects at larger distances 
than this are difficult to constrain because they would likely 
be too faint for KAIT upper limits to be meaningful. For 
example, we searched for pre-SN outbursts in SN 2006gy. 
Although we detected no pre-SN outbursts, the upper limit 
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Figure 11. Same data format as in Fig. [3] but here we display only zoomed data for the four SNe which have limits within 75 days of 
their respective peaks. We compare the data for each SN to those of SN 2009ip and exclude r) Carinae in this comparison figure. The 
data are plotted relative to the observed peak magnitude in published data for each of the SNe if the light curve is well sampled. If our 
KAIT light curve provides a better constraint on the date of the peak (as in the case of SN 2008fq), then we use it instead. 
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of M « —15 mag at d = 73.1 Mpc is not very constraining. 
Ideally, a survey geared toward detecting SN Iln progeni¬ 
tor outbursts should strive to achieve a limiting magnitude 
of m « 21.5 in order to be able to reach sensitivities of 
M ~ —12 mag out to 50 Mpc. The Large Synoptic Survey 
Telescope (LSST) project aims to reach m « 27 mag lim¬ 
its in stacked images (mag 24 in single 15 s exposures) 0 At 
this level of sensitivity, we could measure the progenitors to 
these SNe Iln in their relatively quiescent states if they are 
as luminous as SN 2009ip or rj Carinae out to d « 500 Mpc, 
but variability is needed to pick them out of the confusion- 
limited background. If these objects are not detected by such 
a survey, then it is likely they are being obscured by cir- 
cumstellar dust, in which case an infrared survey would be 
interesting. Even so, deeper optical surveys would allow us 
to statistically compare the number of SNe Iln that arise 
from objects that undergo tremendous mass loss just prior 
to explosion. 
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